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Abstract 

We investigate parton-branching methods based on transverse-momentum dependent (TMD) 
parton distributions and matrix elements for the Monte Carlo simulation of multi-particle final 
states at high-energy colliders. We observe that recently measured angular correlations in ep 
final states with multiple hadronic jets probe QCD coherence effects in the space-like branching, 
associated with finite-angle gluon radiation from partons carrying small longitudinal momenta, and 
not included in standard shower generators. We present Monte Carlo calculations for azimuthal 
two-jet and three-jet distributions, for jet multiplicities and for correlations in the transverse- 
momentum imbalance between the leading jets. We discuss comparisons with current experimental 
multi-jet data, and implications of corrections to collinear-ordered showers for LHC final states. 
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I. INTRODUCTION 



Hadronic final states containing multiple jets have been investigated at the Tevatron and 
HERA colliders, and will play a central role in the Large Hadron Collider (LHC) physics 
program. The interpretation of experimental data for such final states relies both on pertur- 
bative multi-jet calculations (see [1] for a recent overview) and on realistic event simulation 
by parton-shower Monte Carlo generators (see e.g. [2-4]). 

Owing to the complex kinematics involving multiple hard scales and the large phase space 
opening up at very high energies, multi-jet events are potentially sensitive to effects of QCD 
initial-state radiation that depend on the finite transverse-momentum tail of partonic matrix 
elements and distributions. For an overview see [5]. In perturbative multi-jet calculations 
truncated to fixed order in as [1], finite- contributions are taken into account partially, 
order-by-order, through higher-loop corrections. This is generally sufficient for inclusive jet 
cross sections, but likely not for more exclusive final-state observables. 

On the other hand, standard shower Monte Carlos reconstructing exclusive events, such 
as Herwig [6] and Pythia [7], are based on coUinear evolution of the initial-state jet. 
Finite-A;_L contributions are not included, but rather correspond to corrections [8-10] to the 
angular or transverse-momentum ordering implemented in the parton-branching algorithms. 
The theoretical framework to take these corrections into account is based on using initial- 
state distributions (pdfs) unintegrated in both longitudinal and transverse momenta [10], 
coupled to hard matrix elements (ME) suitably defined off mass shell. See e.g. [11] for 
discussion of the Monte Carlo shower implementation of the method. Event generators 
based on /c_|_-dependent showers of this kind include [12-15]. These generators are not as 
developed as standard Monte Carlos like Herwig and Pythia. However, they have the 
potential advantage of a more accurate treatment of the space-like parton shower at high 
energy. 

CoUinear-based shower generators like Herwig and the new Pythia contain the effects 
of color coherence for soft gluon emission from partons carrying longitudinal momentum 
fraction x ~ 0{1). However as the energy increases and emissions that are not coUinearly 
ordered become more important, coherence effects from space-like partons carrying momen- 
tum fractions x <^ 1 set in. These small-x coherence effects are not included in Herwig or 
Pythia but are included in /c^-dependent parton showers, and characterize the structure 
of the initial-state branching at very high energies. 

This paper examines how corrections to space-like parton showers affect properties of 
final state jet correlations and associated distributions. We study azimuthal correlations and 
transverse-momentum correlations for multi-jet processes. We obtain numerical Monte Carlo 
results for coUinear and A;_L-dependent parton showers, and use the precise experimental data 
on tri-jets in ep collisions that have recently become available [16]. We observe significant 
corrections arising from regions [5] with three well-separated hard jets in which the partonic 
lines along the decay chain in the initial state are not ordered in transverse momentum. 
These give rise to quite distinctive features in the jet angular correlations. 

Besides jet final states, the coherence effects from highly off-shell processes discussed in 
this paper affect a variety of different final states at high energy. A significant example 
concerns the associated production of heavy flavors and heavy bosons at the LHC with two 
high-pi jets. We come back to this at the end of the paper in Sec. V. 

The paper is structured as follows. We begin in Sec. II by describing experimental 
results on multi-jet correlations. In Sec. HI we recall basic aspects on the implementation 
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FIG. 1: Dijet azimuthal correlations measured by DO along with the Herwig and Pythia re- 
sults [17]. 

of transverse-momentum dependent pdfs and MEs in parton-branching algorithms. We 
then compute angular correlations in three-jet final states by /c^-dependent Monte Carlo 
showering. We compare the results with Herwig and with experimental data. We consider 
correlations in the azimuthal angle between the two hardest jets, and further analyze the 
distribution of the third jet. We investigate in particular the quantitative effect of the finite 
high-/cx tail in the hard ME. In Sec. IV we present results for jet multiplicity distributions 
and for momentum correlations. In Sec. V we discuss prospects for LHC final states and 
give conclusions. Some details on u-pdf fits and on time-like showering effects are left to 
Appendix A and Appendix B. 

II. MEASUREMENTS OF FINAL-STATE JET CORRELATIONS 

In this section we recall experimental results from Tevatron and HERA on angular cor- 
relations in multi-jet production. 

In a multi-jet event the correlation in the azimuthal angle A0 between the two hardest 
jets provides a useful measurement, sensitive to how well QCD multiple-radiation effects 
are described. In leading order one expects two back-to-back jets; higher-order radiative 
contributions cause the distribution to spread out. At the LHC, measurements of Ac;/) 
distributions in multi-jet events may become accessible relatively early, and be used to test 
the Monte Carlo description of the events. 

Fig. 1 [17] shows the Tevatron A0 measurements. The data are compared with Herwig 
and Pythia results. The data are found [17, 18] to have little sensitivity to final-state 
showering parameters and to be in contrast very sensitive to initial-state showering param- 
eters. In particular, they have been used for re-tuning of these parameters in Pythia [18]. 
A reasonably good description of the measurements by Monte Carlo is obtained. 

On the other hand, the HERA A0 measurements [16, 19, 20] are not so well described 
by the standard Herwig and Pythia Monte Carlo showers in most of the data kinematic 
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range. We will discuss more on this below. These measurements are characterized by the 
large phase space available for jet production and relatively small values of the ratio between 
the jet transverse momenta and center-of-mass energy. For these reasons, despite the much 
lower energy at HERA, they may be just as relevant as the Tevatron data for extrapolation 
of initial-state showering effects to the LHC. 

In the rest of this section we focus on the recent, precise ep measurements [16] of jet 
correlations, and discuss potential sources of large QCD corrections. 

In Ref. [16] the ZEUS collaboration have presented data for two-jet and three-jet produc- 
tion associated with 

Q2 > 10 GeV^ , 10"^ <x< 10-2 , (1) 

and performed a comparison with next-to-leading-order calculations [21]. ZEUS measured 
differential distributions as functions of jet transverse energy and pseudorapidity as well 
as correlations in azimuthal angles and transverse momenta. The selection cuts on the jet 
phase space are given by 

4'HCA/>7GeV , 4^^5'^, >5GeV , -1< rjiab < 2.5 , (2) 

where Et^hcm are the jet transverse energies in the hadronic center-of-mass frame, and rjiah 
are the jet pseudorapidities in the laboratory frame. The overall agreement of data with 
NLO results is within errors [16]. However, while inclusive jet rates are reliably predicted 
by NLO perturbation theory, jet correlations turn out to be affected by large theoretical 
uncertainties at NLO. Results from [16] for di-jet distributions are reproduced in Fig. 2 for 
easier reference. 

The plot at the top in Fig. 2 shows the x-dependence of the di-jet distribution integrated 
over A(/) < 27r/3, where A0 is the azimuthal separation between the two high-E'r jets. The 
plot at the bottom shows the di-jet distribution in for different bins of x. We see that 
the variation of the predictions from order-a^ to order-a^ is significant. In the azimuthal 
correlation for a given x bin, the variation increases with decreasing A0. In the distribution 
integrated over A</), the variation increases with decreasing x. The lowest order, where the 
differential cross section da/d^cf) is non-trivial, is 0{a1) and the NLO calculation is labeled 
with 0{al). 

Given the large difference between order-a^ and order-af results, it seems to be ques- 
tionable to estimate the theoretical uncertainty at NLO from the conventional method of 
varying the renormalization/factorization scale. 

Besides angular distributions, a behavior similar to that described above is also found 
in [16] for other associated distributions such as momentum correlations.^ We will come 
back to this in Sec. IV. 

Note that the Tevatron A0 distribution in Fig. 1 drops by two orders of magnitude 
over a fairly narrow range, essentially still close to the two-jet region. The measurement is 
dominated essentially by leading-order processes. Not surprisingly the Monte Carlos provide 
a good description of the data. In Fig. 2 a comparable two order of magnitude drop occurs 
over the whole A0 range. Much more QCD dynamics than leading order is probed. 

The stability of predictions for the jet observables under consideration in Fig. 2 depends on 
a number of physical effects. Part of these concern the jet reconstruction and hadronization. 
The ZEUS [16] and HI [19, 20] jet algorithm has moderate hadronization corrections [22] 



On the other hand, NLO results are much more stable m the case of inclusive jet cross sections [16]. 
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FIG. 2: (top) Bjorken-x dependence and (bottom) azimuth dependence of di-jet distributions at 
HERA as measured by ZEUS [16]. 



and is free of nonglobal single-logarithmic components [23]. The kinematic cuts [16] on the 
hardest jet transverse momenta are set to be asymmetric, so as to avoid double-logarithmic 
contributions in the minimum pt [24]. Note that > 10 GeV^ (Eq. (1)), and nonper- 
turbative corrections affecting the jet distributions at the level of inverse powers of Q are 
expected to be moderate. 

Further effects concern radiative corrections at higher order. Fixed-order calculations 
beyond NLO are not within present reach for multi-jet processes in ep and pp collisions. Re- 
summed calculations of higher-order logarithmic contributions from multiple infrared emis- 
sions are performed with next-to-leading accuracy in [25] . These contributions are enhanced 
in the region where the two high-i?T jets are nearly back-to-back. Multiple infrared emis- 
sions are also taken into account by parton-branching methods in shower Monte Carlos such 
as Herwig [6]. Note however that important corrections in Fig. 2 arise for decreasing A(^, 
where the two jets are not close to back-to-back and one has effectively three well-separated 
hard jets [5]. The corrections increase as x decreases. Effects analogous to those in Fig. 2 
are seen in the ZEUS results for the three-jet cross section, shown in Fig. 3 [16], particu- 
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FIG. 3: Three-jet cross section versus azimuthal separation between the two highest-E'^^ jets as 
measured by ZEUS [16]. 



larly for the small-A0 and small-x bins.^ In Sec. Ill we analyze the angular distribution of 
the third jet, and find significant contributions for small A0 from regions of the space-like 
shower where the transverse momenta in the initial-state decay chain are not ordered. These 
contributions are not fully taken into account either by fixed-order calculations truncated to 
NLO or by parton showers implementing collinear ordering such as Herwig and Pythia. 

In the next section we present the results of computing jets' angular correlations by 
parton-shower methods that include finite-Zc^ corrections to collinear ordering. We compare 
these results with the collinear-based shower Herwig, and with experimental data. 



III. ANGULAR CORRELATIONS FROM K_l SHOWER MONTE CARLO 



Corrections to the collinear ordering in the space-like shower can be incorporated in 
Monte Carlo event generators by implementing transverse-momentum dependent (TMD) 
parton distributions (unintegrated pdfs) and matrix elements (ME) through high-energy 
factorization [10]. This method allows parton distributions at fixed k± to be defined gauge- 
invariantly for small x. Basic aspects of the parton-shower implementation of the method 
are discussed in [11]. In this section we start by briefly recalling the basis for the introduction 
of unintegrated pdfs (u-pdfs) at high energy; we comment on generalizations relevant for 
low energies and general-purpose tools; then we apply the fc^-dependent parton branching 
to the study of angular jet correlations. 



^ The error band for the theory curves in Fig. 3 [16] is obtained by varying the value of the renormahzation 
scale from (Q^ + Ej,) to (Q^ + Erp)/1Q^ where Et is the average Et of the three hardest jets in each 
event. 
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A. Unintegrated pdfs 

To characterize a transverse-momentum dependent parton distribution gauge-invariantly 
over the whole phase space is a nontrivial question [26, 27], currently at the center of much 
activity. See overview in [5] . In the case of small x a well-prescribed, gauge- invariant defi- 
nition emerges from high-energy factorization [10] , and has been used for studies of collider 
processes both by Monte Carlo (see reviews in [5, 28]) and by semi- analytic resummation 
approaches (see [29, 30]). 

The diagrammatic argument for gauge invariance, given in [10], and developed in [31], 
is based on relating off-shell matrix elements with physical cross sections at x <^ 1, and 
exploits the dominance of single gluon polarization at high energies.^ The main reason why 
a natural definition for TMD pdfs can be constructed in the high-energy limit is that one 
can relate directly (up to perturbative corrections) the cross section for a physical process, 
say, photoproduction of a heavy-quark pair, to an unintegrated gluon distribution, much 
as, in the conventional parton picture, one does for DIS in terms of ordinary (integrated) 
parton distributions. On the other hand, the difficulties in defining a TMD distribution over 
the whole phase space can largely be associated with the fact that it is not obvious how to 
determine one such relation for general kinematics. 

The evolution equations obeyed by TMD distributions defined from the high-energy limit 
are of the type of energy evolution [32] . Factorization formulas in terms of TMD distribu- 
tions [10] have corrections that are down by logarithms of energy rather than powers of 
momentum transfer. On the other hand, it is important to observe that this framework al- 
lows one to describe the ultraviolet region of arbitrarily high k^ and in particular re-obtain 
the structure of QCD logarithmic scaling violations [29-31]. This ultimately justifies the use 
of this approach for jet physics. In particular it is the basis for using corresponding Monte 
Carlo implementations [11-15, 33] to treat multi-scale hard processes at the LHC. 

From both theoretical and phenomenological viewpoints, it is one of the appealing features 
of the high-energy framework for TMD distributions that one can relate its results to a well- 
defined summation of higher-order radiative corrections. By expanding these results to fixed 
order in a^, one can match the predictions thus obtained against perturbative calculations. 
This has been verified for a number of specific processes at next-to-leading order (see for 
instance [34] for heavy fiavor production) and more recently at next-to-next-to-leading order 
(see for instance [35]). Note that this fact also provides the basis for shower algorithms 
implementing this framework to be combined with fixed-order NLO calculations by using 
existing techniques for such matching. 

Later in this section we use Monte Carlo implementing the high-energy definition of u- 
pdfs to analyze jet production. Before doing this, we comment briefiy on open issues and 
generalizations to low energies. 

B. Comments on unintegrated pdfs beyond low x 

In the general case, factorization formulas in terms of unintegrated parton distributions 
will have a considerably complex structure [26]. Full results are yet to be established. A 



^ It is emphasized e.g. in [27, 28] that a fully worked out operator argument, on the other hand, is highly 
desirable but is still missing. 
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prototypical calculation that illustrates this structure is carried out in [36], which treats, 
rather than a general scattering observable, a simpler problem, the electromagnetic form 
factor of a quark. This case is however sufficient to illustrate certain main features, in 
particular the role of nonperturbative, gauge-invariantly defined factors associated with 
infrared subgraphs (both collinear and soft), and the role of infrared subtractive techniques 
that serve to identify these factors. See also [37] for recent analyses along these lines for 
more general processes involving fully unintegrated pdfs. 

One of the questions that a full factorization statement will address is the treatment of 
soft gluons exchanged between subgraphs in different collinear directions. The underlying 
dynamics is that of non-abelian Coulomb phase, treated a long time ago in [38] for the 
fully inclusive Drell-Yan case. But a systematic treatment for more complex observables, 
including color in both initial and final states, is still missing, as emphasized recently in [39- 
41] for di-hadron and di-jet hadropro duct ion near the back-to-back region.^ 

A further question concerns lightcone divergences [26] and the x ^ 1 endpoint behavior. 
The singularity structure at x — 1 is different in the TMD case than for ordinary (integrated) 
distributions, giving divergences even in dimensional regularization with an infrared cut- 
off [43] . The singularities can be understood in terms of gauge-invariant eikonal-line matrix 
elements [43], and the TMD behavior can be related to cusp anomalous dimensions [44, 45] 
and lack of complete KLN cancellations [46]. In general this affects the precise form of 
factorization and relation with collinear distributions. 

Applications of u-pdfs at low energies include semi-inclusive leptoproduction ([47-49], and 
references therein), spin asymmetries [50], exclusive reactions [51]. In these cases infrared 
subtractive techniques of the type [36, 52] serve for TMD-factorization calculations [53] and 
in particular for the proper treatment of overlapping momentum regions.^ At high-energy 
colliders, general characterizations of TMD distributions will be relevant for turning present 
k^-showering generators into general-purpose tools to describe hadronic final states over the 
whole phase space [5, 60]. 

In the rest of this section we will consider applications of k^-shower generators to multi- 
jet final states. The main focus is on regions where jets are far from back-to-back, and 
the total energy is much larger than the transferred momenta so that the values of x are 
small. In this regime the ambiguities related to soft Coulomb exchange and to lightcone 
divergences are not expected to be crucial. We will find that the TMD distributions, as 
well as the transverse-momentum dependence of short-distance matrix elements, play a very 
essential role to describe correlations in angle and momentum of the jets. 

C. k± shower with u-pdfs 

Monte-Carlo event generators based on unintegrated pdfs use factorization at fixed k_|_ [10] 
in order to a) generate the hard scattering event, including dependence on the initial trans- 



^ Note that interestingly in [42], which has a different point of view than TMD, Coulomb/radiative mixing 
terms are found to be responsible for the breaking of angular ordering in the initial-state cascade and the 
appearance of superlcading logarithms in di-jet cross sections with a gap in rapidity. 

^ Subtraction techniques related to those of [36, 52] are developed in [54] for soft-coUincar effective theory, 
and studied in [55] and [56] in relation with standard perturbativc methods. Sec also SCET applications 
to shower algorithms [57], TMD pdfs [58] and jet event shapes [59] for use of these techniques. 
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verse momentum, and b) couple this to the evolution of the initial state to simulate the 
parton cascade. Implementations of this kind include [8, 12-15]. The hard scattering event 
is generated by A;_L-dependent matrix elements (ME) computed from perturbation theory. 
Different generators differ by the detailed model for initial state. For the calculations that 
follow we use the Monte Carlo implementation Cascade [12]. 

The hard ME in the Monte Carlo are obtained by perturbative calculation [10], while the 
u-pdfs are determined from fits to experimental data [33]. The parton-branching equation 
used for the unintegrated gluon distribution A is schematically of the form [8, 12, 33] 

A{x,k±,n) = Ao{x,k±, n) + J J -^'^(i^-^q) 

X A{fi,zq)P{z,q,k^) Ai^,k^ + {1- z)q,q) . (3) 

The first term in the right hand side of Eq. (3) is the contribution of the non-resolvable 
branchings between starting scale Qq and evolution scale /i, and is given by 

Ao{x,k±,n) = Ao{x,k:i,Qo) A{n,Qo) , (4) 

where A is the Sudakov form factor, and the starting distribution Ao{x, k±, Qq) at scale Qq 
is determined from data fits. Details on the starting distribution used for the calculations 
that follow are given in Appendix A. 

The integral term in the right hand side of Eq. (3) gives the /c^-dependent branchings 
in terms of the Sudakov form factor A and unintegrated splitting function V. The explicit 
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FIG. 4: (left) Coherent radiation in the space- like parton shower for x <^ 1; (right) the unintegrated 
splitting function V, including small-x virtual corrections. 

expressions for these factors are specified in [33], and include the effects of coherent gluon 
radiation not only at large x (as e.g. in Herwig) but also at small x [9] in the angular 
region (Fig. 4) 

a/x > ai > a , (5) 

where the angles a for the partons radiated from the initial-state shower are taken with 
respect to the initial beam jet direction, and increase with increasing off-shellness. Unlike 
conventional showers, the splitting function V depends on transverse momenta and includes 
part of the virtual corrections, in such a way as to avoid double counting with the Sudakov 
form factor while reconstructing color coherence in the small-x region (5). 

The Monte Carlo also contains time-like parton showering. The impact of time-like 
showers on the jet observables that we will examine in this section turns out to be very 
small. This is similar to what is found in the studies [17, 18] of Tevatron di-jets and A0 
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distribution, based on Herwig and Pythia. See remark toward the beginning of Sec. II. 
Some details on the treatment of time-hke showering effects are reported in Appendix B. A 
more complete account of this topic may be found in [61]. 



D. Azimuthal jet distributions 

The /c_L-dependent ME and parton branching lead to a different angular pattern of initial- 
state gluon radiation than standard, collinear-based showers, e.g. Herwig. In particular, 
while the Herwig angular ordering reduces to ordering in transverse momenta for x — > 0, 
the /c^-dependent shower contains finite-angle corrections in this limit [11]. We now com- 
pute angular distributions for the ep three-jet cross section by the /c^-shower Monte Carlo 
Cascade and by Herwig. Let A0 be the azimuthal separation between the two jets with 
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FIG. 5: Cross section in the azimuthal angle A(j)i3 between the hardest and the 3rd jet for small 
(A(/> < 2, top) and large (Ac/) > 2, bottom) azimuthal separations between the leading jets. The 
A;_LMonte Carlo results Cascade are compared with Herwig. 
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the highest transverse energy Et-, 



yjct-l 



fjct-2 , 



(6) 



where the azimuthal angle for each jet is defined in the hadronic center-of-mass frame. 
Similarly, we define Ac^is as the azimuthal separation between the hardest and the third jet. 

In Fig. 5 we compute the three-jet cross section and plot it versus the azimuthal angle 
A013, by distinguishing the cases in which the two leading jets are at small angular separation 
(A0 < 2) or large angular separation (Ac/) > 2). Cascade gives large differences from 
Herwig in the region where the azimuthal separations A0 between the leading jets are 
small, see top plot of Fig. 5. This refiects the fact that at small A0 the phase space opens 
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are multiplied by a factor of 2. 



up for events in which the partonic lines along the initial decay chain are not ordered in 
transverse momentum. Such configurations are taken into account in Cascade with the 
appropriate matrix element, at least for small enough x, but not in Herwig. The x values 
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considered in Fig. 5 are those corresponding to the three-jet measurements in [16]. As 
increases, the results from CASCADE and Herwig become closer. See bottom plot of Fig. 5. 
This is associated with the fact that for A0 approaching the back-to-back region the phase 
space for finite-/cx emissions is reduced. In this region one thus expects both Monte Carlos 
to give reasonable approximations. 

Fig. 6 shows the angular correlations for final states with two jets and three jets. We 
compute the azimuthal distribution of di-jet and three-jet cross sections in the separation A0 
between the leading jets. We show the distributions obtained by Cascade and by Herwig, 
compared with the measurement [16]. We multiply the Herwig result by a constant factor 
equal to 2, which the top plot in Fig. 6 shows is the K- factor needed in order to get the 
normalization approximately correct for the two-jet region. Observe that the shape of the 
distribution is different for the two Monte Carlos. As expected from the result of Fig. 5, 
Cascade gives the largest differences to Herwig at small A0, and becomes closer to 
Herwig as A0 increases. In particular, we observe that while the K-factor of 2 for Herwig 
is sufficient for the two-jet region, the shape of the jet distribution is not properly described 
by Herwig as A0 decreases. The description of the measurement by Cascade is good, 
whereas Herwig is not sufficient to describe the measurement in the small A0 region. We 
further see in the bottom plot of Fig. 6 that the three-jet cross section is reasonably well 
described by the /c_L-shower result but not by Herwig. 
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Note that the interpretation of the angular correlation data in terms of corrections to 
coUinear ordering is consistent with the finding [16] discussed in Sec. II that while inclusive 
jet rates are reliably predicted by NLO fixed-order results, NLO predictions are affected 
by large corrections to di-jet azimuthal distributions (going from ^(a^) to 0{al)) in the 
small-A(;/) and small-x region, and begin to fall below the data for three-jet distributions in 
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the smallest bins (Fig. 3 [16]). 

The physical picture underlying the /c^-shower calculation in Figs. 5,6 involves both 
transverse-momentum dependent parton distributions (determined from experiment) and 
matrix elements (computed perturbatively). Fig. 7 illustrates the relative contribution of 
these different components to the result, showing different approximations to the azimuthal 
dijet distribution normalized to the back-to-back cross section. The solid red curve is the 
full result. The dashed blue curve is obtained from the same unintegrated pdf's but by 
taking the coUinear approximation in the hard matrix element, 

M{k^) ^ M,omn{kA_) = A^(Ox) e(/i - k^) . (7) 

The dashed curve drops much faster than the full result as A0 decreases, indicating that 
the high-/c^component in the hard ME [10] is necessary to describe jet correlations for small 
A0 [62]. For reference we also plot, with the dotted (violet) curve, the result obtained from 
the unintegrated pdf without any resolved branching, 

A{X, k^, /i) Ano-res.{x, k^, /i) = Ao{x, k±, Qq) A(/i, Qq) . (8) 

Here is the starting distribution at Qo and A is the Sudakov form factor, giving the 
no-radiation probability between Qq and /x. This represents the contribution of the intrinsic 
/cx distribution only, corresponding to nonperturbative, predominantly low-Zc^modes. That 
is, in the dotted (violet) curve one retains an intrinsic k±y^ but no effects of coherence. 
We see that the resulting jet correlations in this case are down by an order of magnitude. 

The results of Fig. 7 illustrate that the /c_L-dependence in the unintegrated pdf alone is 
not sufficient to describe jet production quantitatively, and that jet correlations are sensitive 
to the finite, high-Zc^tail of matrix elements [10] computed from perturbation theory. We 
note that the inclusion of the perturbatively computed high-k^ correction distinguishes the 
present calculation of multi-jet cross sections from other shower approaches (see e.g. [15]) 
that include transverse momentum dependence in the pdfs but not in the matrix elements. 
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FIG. 8: Comparison of the k_L-shower Cascade with the NLO di-jet calculation Disent: (left) dis- 
tribution in single-jet transverse energy; (right) distribution in the di-jet transverse energy. 



To examine more closely the distribution in k^ that results from highly off-shell sub- 
processes, in Fig. 8 we study the jet cross section in transverse energy and compare the 
k^-shower with the NLO result from the Disent event generator. It is noteworthy that 
the large-p( part of the di-jet spectrum is very close for the two calculations. At low pt one 
sees the Sudakov form-factor effect in the shower result. Differences in the single-jet spectra 
are also of interest and can be shown to be associated with quark contributions [61]. These 
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detailed comparisons may be of use to relate [63] DIS event shapes measuring the transverse 
momentum in the current region to hadro-production px spectra. 
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FIG. 9: Azimuthal di-jet distribution obtained from the expansion of the k_L-shower Cascade to 
one-ghion emission level, compared with the NLO di-jet calculation DiSENT(gluon channel). 



In Fig. 9 we push further the comparison at next-to-leading-order level. We switch off 
hadronization, and use the k^-shower Monte Carlo Cascade as a parton-level generator. 
We evaluate the k^-shower by expanding in the number of extra emissions, and truncate to 
the level of one gluon emission. We compare this with the NLO Disent calculation, taking 
only the gluon channel in DiSENT. We compute the azimuthal di-jet distribution at various 
values of x. The plots in Fig. 9 indicate that for sufficiently small x the one-gluon expansion of 
the shower program agrees with the full NLO result. We view this as a numerical consistency 
check of the shower program in the case of a relatively complicated final-state correlation, 
to be considered jointly with the analytic cross-checks quoted in Sec. Ill A, e.g. [34], for the 
case of analytic small-x results for inclusive observables. 

We conclude this section by observing that the jets that we are considering are produced 
in the region of rapidities of Eq. (2), away from the forward region. While forward- region 
observables are relevant in their own right and have long been studied as probes of the 
initial-state shower dynamics (see e.g. [4, 28] and references therein), Monte Carlo results 
for such observables have a more pronounced dependence on the details of the model used 
for u-pdf evolution [11] (see also discussions in [27, 37, 48]). It is thus interesting that 
significant effects of non-ordering in k± for the space-like shower are found in the present 
case for centrally produced hard jets. 



IV. JET MULTIPLICITIES AND MOMENTUM CORRELATIONS 

We now turn to jet multiplicities and transverse-momentum correlations. These observ- 
ables provide further details on the structure of the multi-jet final states. As noted in Sec. II, 
several of the transverse-momentum correlations measured in [16] are affected by sizeable 
theory uncertainties at NLO [16, 21]. 

Let us first consider jet multiplicity distributions. Finite-Zc^ corrections increase the mean 
gluon multiplicity and broaden the spectrum [8-10]. In Fig. 10 we compute the distribution 
in the number of jets A^, normalized to the two-jet cross section cr. As in Fig. 5 we show 
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FIG. 10: Jet multiplicities obtained by Cascade and Herwig for (top) A(j) < 2 and (bottom) 
A<p > 2. 



separately the results for small and large azimuthal separations between the hardest jets, 
A0 < 2 and A0 > 2. Jet multiplicities at small Acf) are where the clearest differences 
appear between the two parton showers. The fc^-shower result receives larger contribution 
from high multiplicities. Besides the absolute size of this contribution, note that Fig. 10 
illustrates the difference in the shape between the two Monte Carlos. 

In Ref. [16] the ZEUS collaboration has presented measurements of various momentum 
correlations. We examine two such distributions for three-jet cross sections in Figs. 11 and 
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FIG. 11: Momentum correlations obtained by Cascade and Herwig, compared with ep data [16]: 
three-jet cross section versus the variable | introduced in the text. 

12. In Fig. 11 is shown the distribution in the magnitude of the sum of the transverse 
momenta pt for the two jets with the highest E^, \ YIPt'^I- The back-to-back region corre- 
sponds to I Y2Pt^\ — ^ in this plot. The region of large | XIPt^I is the region with at least 
three well-separated hard jets. The /c_L-shower result describes this region reasonably well. 
The results from Herwig are quite lower. 

In Fig. 12 is the distribution in the vector difference of the highest-i?T jet transverse 
momenta, scaled by twice the transverse energy of the hardest jet, \Ap}j^'^\/{2E^). The back- 
to-back region corresponds to |Ap^ \/{2E^) — > 1 in this plot. The behavior of the Monte 
Carlo results compared to the data is rather similar to that in Fig. 11. 

In summary, the calculations of this paper show that the fc^-shower results describe 
well the shape of multi-jet distributions observed experimentally, including correlations, and 
give quite distinctive features of the associated distributions compared to standard showers 
such as Herwig. The largest differences between the two parton showers occur when the 
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FIG. 12: Momentum correlations obtained by Cascade and Herwig, compared with ep data [16]: 
three-jet cross section versus the variable | Ap;^^|/(2ii^-^) introduced in the text. 

azimuthal separations between the leading jets are small, whereas the results become more 
similar in the two-jet region. See e.g. Figs. 6,10. In the region of small azimuthal distances 
the largest variation occurs between order-a^ and order-a^ results in the fixed-order NLO 
calculations, particularly for small x. In cases where corrections are not large, the NLO 
and /c_L-shower calculations are rather close. The results support a physical picture of multi- 
jet correlations in which sizeable radiative corrections arise not only from collinear/soft 
emission, included in Herwig as well, but also from finite-angle emission, associated with 
the growth of transverse momenta transmitted along the space-like jet. Small-x coherence 
effects, computed in this section and the previous section for jet multiplicities, momentum 
correlations and angular correlations, are included in the /c^-shower but not in Herwig. 
They are associated with multi-gluon radiation terms to the initial-state shower that become 
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non- negligible at high energy and small A0.^ 

The above observations suggest the usefulness of combining NLO and /c_L-shower for a 
broad range of multi-jet observables, in order to obtain more reliable predictions over a wider 
kinematic region. Monte Carlo results depend on the maximum angle parameter fi [8, 11, 12] 
at which the shower is evaluated. The perturbative matching will involve this angle. Studies 
of the dependence of Monte Carlo results on /_i will be reported elsewhere. 

V. PROSPECTS FOR LHC FINAL STATES AND CONCLUSIONS 

Experimental analyses of multi-particle final states at the Large Hadron Collider depend 
on realistic parton-shower Monte Carlo simulations. Multi-particle production acquires qual- 
itatively new features at the LHC compared to previous hadron-hadron experiments due to 
the large phase space opening up for events characterized by multiple hard scales, possibly 
widely disparate from each other. This brings in both potentially large radiative corrections 
and potentially new effects in the nonperturbative components of the process being probed 
near phase-space boundaries. It is not at all obvious that the approximations involved in 
standard Monte Carlo generators that have successfully served for event simulation in past 
collider experiments will be up to the new situation. 

In this paper we have discussed the method of k^-dependent Monte Carlo shower, based 
on transverse-momentum dependent (TMD), or unintegrated, parton distributions and ma- 
trix elements defined by high-energy factorization. The main advantage of the method over 
standard Monte Carlo generators is the inclusion of corrections to collinear-ordered showers, 
and of effects of QCD coherence associated with finite-angle radiation from space-like par- 
tons carrying arbitrarily soft longitudinal momenta. Sensitivity to these dynamical features 
is bound to be enhanced by the high-energy multi-scale kinematics. The theoretical basis 
of the k^-shower method allows one to go to arbitrarily high transferred-momentum scales, 
thus making it suitable for event simulation of jet physics at the LHC. 

In the paper we have pointed to developments of the approach toward general-purpose 
event generators, and illustrated validation of k^-shower Monte Carlo using experimental ep 
data for final states with multiple hadronic jets. We have noted that while Tevatron di-jet 
correlations are dominated by leading-order processes, and are reasonably well described by 
coUinear-based event generators, this is not so in the case of ep data. We have found that 
including finite-k^ radiative contributions in the initial state shower gives sizeable effects 
and improves significantly the description of angular correlations and transverse-momentum 
correlations. Despite the lower ep energy, the multi-jet kinematic region considered is char- 
acterized by the large phase space available for jet production and relatively small values of 
the ratio between the jet transverse momenta and center-of-mass energy, and is thus relevant 
for extrapolation of initial-state showering effects to the LHC. 

Besides jet final states, the corrections to collinear-ordered showers that we are treating 
will also affect heavy mass production at the LHC, including final states with heavy bosons 
and heavy flavor. An example is provided by bottom-quark pair production. Going from 
the Tevatron to the LHC [64] implies a sharp increase in the relative fraction of events domi- 
nated by the g ^ bb subprocess coupling [10] to the spacelike jet. This is bound to affect the 



^ Near the back-to-back region of large A(j>, on the other hand, corrections due to mixed Coulomb/radiative 
terms can also become important [39, 40]. 
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reliability of shower calculations based on collinear ordering (as well as the stability of NLO 
perturbative predictions), as these do not properly account for contributions of bb in associa- 
tion with two hard jets, with pt of the heavy quark pair large compared to the bottom-quark 
mass but small compared to the transverse momenta of the individual associated jets. These 
kinematic regions are the analogue of the regions unordered in k^ studied in this paper for 
jet correlations. The fraction of bb events of this kind is not very significant at the Tevatron 
but will be sizeable at the LHC. The quantitative importance of unordered configurations 
coupling to g —>■ bb will reduce the numerical stability of coUinear-based predictions (NLO, 
or parton-shower, or their combination [65]) with respect to renormalization/factorization 
scale variation in the case of LHC. On the other hand, these are precisely the configurations 
that the k^ Monte Carlo shower is designed to treat. 

Even more complex multi-scale effects are to be expected, and are beginning to be in- 
vestigated [61], in the associated production of bottom quark pairs and W/Z bosons [66], 
and possibly in final states with Higgs bosons [33, 67]'^ especially for measurements of the 
less inclusive distributions and correlations. The vector boson case is relevant for early phe- 
nomenology at the LHC, as small- a; broadening of W and Z px distributions [69] (see [70]) 
affects the use of these processes as luminosity monitor [71]. 

The kx-shower method discussed in this paper can be used all the way up to high 
transferred-momentum scales. As an illustration in Fig. 13 we present a numerical cal- 
culation for the transverse momentum spectrum of top-antitop pair production at the LHC. 
Small-x effects are not large in this case. Rather, this process illustrates how the shower 
works in the region of finite x and large virtualities on the order of the top quark mass. 
It is interesting to note that even at LHC energies the transverse momentum distribution 
of top quark pairs calculated from the k^-shower is similar to what is obtained from a full 
NLO calculation (including parton showers, MC@NLO [65]), with the k^-shower giving a 
somewhat harder spectrum. Fig. 13. 
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FIG. 13: Comparison of transverse momentum distribution of tt pairs calculated from the k^- 
shower CASCADE with the NLO calculation MC@NLO at LHC energies. 

We conclude by observing that using off-shell matrix elements convoluted with unin- 



Non-negligible numerical effects of high-energy subleading terms were noted [68] in the predictions for the 
Higgs transverse- momentum spectrum at the LHC. 
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tegrated parton distributions including explicit parton showering, many of the subleading 
effects are properly simulated both in ep collisions and at the LHC. We have found that 
multi-jet predictions provide comparable results to NLO calculations, where applicable, and 
are much closer to the measurements in a region where significant higher order contributions 
are expected. The results provide a strong motivation for systematic studies of k_|_-dependent 
parton branching methods. 
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APPENDIX A: FITS TO THE STARTING PDFS 

The branching equations (3), (4) contain the starting gluon distribution Ao{x, k±,QQ) at 
scale Qq. This is determined from fits to experimental data. In this appendix we report 
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FIG. 14: (top) x-dependence and (bottom) fc^-dependence of the unintegrated gluon distribution 
at different values of the evolution scale /i. 

results for this distribution. 

The starting Ao{x, k±, Qq) at scale Qo is parameterized as [11, 33] 

xAoix, k±, Qq) = A x-^ (1 - xf exp [-{k± - . (Al) 
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The values of the parameters A, B, C, A and u in Eq. (Al) are determined from data 
fits [33, 72]. In the calculations of the present paper we use the u-pdf set specified by the 
following parameter values: 

Qo = 1.1 GeV , A = 0.4695 , B = 0.025 , 
C = 4.0 , A = 1.5 GeV , u = (1.5/^2) GeV . (A2) 

In Fig. 14 we plot the x-dependence and /c_L-dependence of the resulting gluon distribution 
at different values of the evolution scale 



APPENDIX B: TIME-LIKE SHOWERING EFFECTS 

The partons from the initial state cascade are allowed to develop a time-like shower in 
Cascade 2.0.2, to be published in [61]. Full details will be reported in this publication. To 
give an idea of the effects, we include one of the results in this appendix. 
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FIG. 15: The effect of including time-like showering on the ep charged-particle spectrum, along 
with the data [73]. 



The maximum scale for the time-like cascade is given by the transverse momentum of 
the initial state gluon. No additional constraints are applied to the time-like shower. It 
is found that the number of gluons after the initial state cascade with time-like showering 
increases; however the effect on the angular correlations considered in Sec. Ill of this paper 
is negligible (and smaller than the statistical error of the Monte Carlo simulation). 
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In observables which are more sensitive to the time-hke showering, hke the charged- 
particle transverse momentum spectra (Fig. 15), a small effect coming from the time-like 
showering can be observed, and is of the same size as that obtained from Monte Carlo event 
generators using the coUinear parton-showering approach. 
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